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The corresponding Hamiltonian is defined as the expectation value of the electromag- 
netic term in the Einstein-Maxwell Hamiltonian constraint, regularized a la Thiemann, 
with respect to a would-be semiclassical state. The resulting energy dispersion relations 
entail Planck scale corrections to those in flat spacetime. Both the helicity dependent 
contribution of Gambini and Pullin [^J and, for a value of a parameter of our approxi- 
mation, that of Ellis et. al. |]lO[ are recovered. The electric/magnetic asymmetry in the 
regularization procedure yields nonlinearities only in the magnetic sector which are briefly 
discussed. Observations of cosmological Gamma Ray Bursts might eventually lead to the 
needed accuracy to study some of these quantum gravity effects. 
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1 Introduction 



It has been recently suggested that quantum gravity effects might be indeed observable |T]|-0]- 
Among the diverse phenomena where these observations could take place we can distinguish 
(see e.g. (i) strain noise induced in gravitational wave detectors, (ii) neutral Kaon 

systems and (iii) energy dependent time of arrival of photons/neutrinos from distant sources. 
Experimental sensitivities in each of the above situations can be argued to be at the edge of 
what is required to reveal quantum gravity phenomena, namely Planck length sensitivities. In 
this work we focus on the third possibility. The idea is to look for modified dispersion relations 
of photons with energy E and momentum k, of the form 



where ^ is a numerical factor of order one and Eqg is an energy scale of order < 10^^ GeV, 
which signals the need of considering the quantum character of gravity. The above expression 
leads to the following modification of the speed of light in vacuum 



which implies a retardation time 



with respect to a signal propagating with speed c. When we consider cosmological distances 
L ^ 10^° ly and an energy scale Eqg ^ 10^^ GeV in , the corresponding values are At ^ 
10-^ s {E ^ 20MeV) and At ^ 10"^ s (E ^ 0.20MeV). In order to detect such effects, an 
experimental time resolution 6t at least of order At is required. Thus, short and intense bursts 
traveling large distances would be the best candidates. Recent observations pointing towards 
the possibility of attaining such conditions are: (i) some Gamma Ray Bursts (GRB) originate 
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at cosmological distances, (~ 10^'' ly) and (ii) sensitivities 5t up to submillisecond scale 
have been achieved in GRB observations |0] and they are expected to improve in future spatial 
experiments 

It is thus timely to investigate whether candidate quantum gravity theories can account for 
([l|). Modifications to Maxwell's equations in vacuum, induced by quantum gravity effects, have 
been calculated by Gambini and Pullin . By considering a semiclassical regime in which the 
electromagnetic field is a classical object whereas space is described by loop quantum gravity, 
they obtained the dispersion relations 

uj±{k) = k{lT'2ilpk) , k=\kl (4) 

where ± labels the helicity of the photon and ^p \{]^^'^ cm is the Planck length. This 
modification, being helicity dependent, yields parity violation and birefringent effects. 



On the other hand, a string theory approach proposed by Ellis et al. [T^ suggests that 
a D-brane recoil in the quantum-gravitational foam induces a distortion in the surrounding 
space, which modifies the photon propagation properties. For a review of this approach see for 



example Ref. [|TT|. The dispersion relations obtained via this second procedure are 



with ^ > 0. They arise from parity conserving corrections to Maxwell's equations which lead 
to a first order ( in l/M/j) helicity independent effect in the dispersion relations that is linear 
in the photon energy. No birefringent effects appear at this order. In this approach, the 
red-shifted difference in the time arrival of two photons with present-day energies Ei and E2 
has been calculated. For the BATSE data, when the redshifts z of the GRB are known, a 
small subset of coincident photon pulses corresponding to channel 1 (20-50 keV) and channel 
3 (100-300 keV) are fitted and At is calculated [|1^. No significant effect in the data available 



is found. On the other hand, none of the pulses studied exhibited microbursts having a short 
time structure on scales < 10^^ s. Were this the case, the sensitivity of the analysis would be 
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greatly improved. Alternative studies based upon effective perturbative quantum gravity 
open system techniques |TB] and quantum light cone ffuctuations |n] have been performed. 

The study of cosmological neutrinos could also provide an excellent arena to probe quantum 
gravity induced propagation effects because space is practically transparent to them, even at 
very high energies. In fact, the fireball model, which is one of the most popular models of GRB, 
predicts the generation of 10^"^ — 10^^ eV Neutrino Bursts (NB) |jT5|, The planned Neutrino 
Burster Experiment (NuBE) will measure the flux of ultra high energy neutrinos (> 10 TeV) 



over a ~ 1km effective area, in coincidence with satellite measured GRB's |Tj]. It is expected 
to detect ~ 20 events per year, according to the fireball model. Other relevant experiment 
aimed at observing ultra high energy cosmic rays, including neutrinos, is the OWL/Airwatch 
project which expects to see ~ 3 x 10^ — 10^ cosmic ray events with energies > lO^'^eV [|T8| , p!9[| . 
Notably, this experiment is able to investigate time correlations among high energy neutrinos 
and gamma-rays. Hence in the foreseeable future, it might be possible to study quantum gravity 
effects on observed astrophysical neutrinos and photons. At the least, such observations could 
be used to restrict quantum gravity theories. 

Motivated by these interesting possibilities, we have calculated the quantum gravity induced 
modifications to neutrino propagation in [^, within the loop quantum gravity framework. We 
obtained corrections to the velocity of propagation which are proportional to {kip) together 
with helicity dependent corrections of order {kipY. The energy dependence in the first case 



coincides with that found later by Ellis et. al. using string theory methods [pT |. 

In this work we extend our approach to the case of photons. The corrections obtained within 
our approximation contain those of [0] and, for a given value of a parameter in our scheme. 



those of |T0|. Besides, we briefly discuss higher order nonlinearities arising in the magnetic 
sector of the effective Maxwell Hamiltonian. 

The organization of the paper is as follows. In section ^ we recall some basic aspects 
of loop quantum gravity which are necessary for our analysis. After reviewing Thiemann's 
regularization of the Hamiltonian constraint of the Einstein-Maxwell theory in section H, we 
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provide a general description of our approximation in section ^. The corrections arising from 
the electric and magnetic sectors are calculated in section ^ Once the effective Hamiltonian is 
obtained, we derive the modified Maxwell equations together with the corresponding dispersion 
relations in section ^. Section |^ contains a brief analysis of the non-linear effects arising from 
the magnetic sector. An outline of red shift effects on cosmological photon time delays is given 
in section |. Finally, section ^ contains a discussion of the results. 

2 Loop quantum gravity 

In this section we summarize the main ingredients defining this approach also denominated 



quantum geometry p2[. Among the main results along this approach one finds: i) well defined 



geometric operators possessing a discrete spectrum, thus evidencing discreteness of space 



ii) a microscopic account for black hole entropy |^ and, more recently, hints on quantum 
avoidance of a would be classical cosmological singularity |2^. (For a review on these topics 



see for example Ref. [pSj] .) 

To begin with it is assumed that the spacetime manifold M has topology x M, with E 
a Riemannian 3-manifold. Here a co-triad is defined, with a,b,c, . . . being spatial tensor 
indices and i,j,k,... being su{2) indices. Thus the corresponding three-metric is given by 
Qab = e^e^. In addition, a field -K"* is defined by Kab = sgn[det{el)]K^g^el, which is related to the 
extrinsic curvature Kab of S . A canonical pair for the gravitational phase space is (i^*, Ej/k), 
where E!^ = ^e'^^^eijkele^ and k is Newton's constant. It turns out that such a canonical pair 
yields a complicated form for the Hamiltonian constraint of general relativity. A convenient 
canonical pair, making this constraint polynomial, was introduced by Ashtekar [Q. Neverthe- 



less, two severe difficulties to proceed with the quantization remained: (i) the implementation 
of a diffeomorphism covariant regularization for the density-weight two Hamiltonian constraint 
hereby obtained and (ii) the extension to non-compact groups of the diffeomorphism covariant 
techniques already developed for gauge theories with compact groups p^. In fact, the Ashtekar 
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variables (*y4Jj = Tl^ — iK^, iEI/n) , with being the torsion free connection compatible 
with e^, are complex valued. Namely the gauge group is SL{2,(P), which is non compact. 

Some proposals to come to terms with difficulty (ii) were: to consider real connection 
variables ||28|, to implement a Wick transform and to define tractable reality constraints 
30| . All of these left open (i). Thiemann subsequently proposed to solve (i) and (ii) by 



incorporating real connection variables while keeping the density weight one character of the 
Hamiltonian constraint. He further provided a quantum version of the theory in the pure 



gravity case, as well as in those cases including the coupling of matter to gravity His 
approach is next reviewed, since we rely upon it for our analysis of the electromagnetic case. 
Let us start with the following canonical pairs for the Einstein-Maxwell system: (A^ = 
+ KI^,E'^/k) for the gravity sector and {AaiE^/Q"^) the electromagnetic sector. The 
latter has gauge group U{1) and Q is the electromagnetic coupling constant, related to the 
dimensionless fine structure constant by = Q^^. The corresponding contributions to the 
Hamiltonian constraint are 

i^Einstcin[iV] = I d^'x N ^=ti {2[K,, K,] - F,,) [E,, E,]) , 
JT. Kvdet q 

i^Maxwell[iV] = I X N —^^^[E^ E' + B\ (6) 

Here Fab is the curvature of Aa and is the magnetic field of the U{\) connection A. The 
actual classical configuration space is the space A/Q ol (both) connections modulo their gauge 
transformations. Indeed, this is what occurs in gauge theories where the fundamental field is a 
connection. This completes the classical description of the phase space of the theory. 

The quantum arena is given as follows As in any quantum field theory, because of 

the infinite number of degrees of freedom, an enlargement of the classical configuration space 
is required. This is far from trivial since the measures defining the scalar product, which are 
required to provide a Hilbert space, get concentrated on distributional fields and hence lie 
outside the classical configuration space. The key idea to build up such an enlargement is to 
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make Wilson loop variables (traces of parallel transport operators) well defined. The resulting 
space A/Q can be thought of as the limit of configuration spaces of lattice gauge theories for 
all possible floating (i.e. not necessarily rectangular) lattices. Hence, geometric structures on 
lattice configuration space are used to implement a geometric structure on A/G- This enables 
to define a background independent calculus on it which, in turn, leads to the construction of 
the relevant measures, the Hilbert space and the regulated operators. 

In line with the Dirac procedure for constrained systems, one first ignores the constraints 
and constructs an auxiliary Hilbert space Tiauxj so that the set of elementary real functions on 
the full phase space is represented by self-adjoint operators in TYaux- It turns out that Tiaux 
is just L^{A/Q, Ho), with iiq being a suitable measure that implements the self-adjointness 
property. 

Diffeomorphism constraints are well defined operators on Tiaux yielding no anomalies. There 
exists a dense subspace $ of TYaux so that its topological dual includes a complete set of 
solutions to the diffeomorphism constraint. They are characterized by generalized knots (i.e. 
diffeomorphism invariant classes of graphs). Besides, a diffeomorphism invariant Hilbert space 
is obtained for such states with an inner product that represents real observables by self-adjoint 
operators. 

Furthermore, Tiaux admits a basis in terms of the so called spin network states. A spin 
network is a triple {a,j,p) consisting of a graph a, a coloring defined by a set of irreducible 
representations (ji, . . . , j„) of SU{2), with ji corresponding to the edge of a and a set of 
contractors (pi, . . . ,Pm)- Here, a contractor is just an intertwining operator from the ten- 
sor product of representations of the incoming edges at the vertex to the tensor product 
of representations of the outgoing edges. Compactness of SU{2) makes the vector space of 
all possible contractors Pk finite, for a given j and vertex Vk. An additional non-degeneracy 
condition is included: jg is not trivial for any edge e and a is taken to be minimal (i.e. any 
other a', occupying the same points in E as a, can always be built by subdividing the edges of 
a, but not the other way around). 
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A spin network state is a cylindrical function (a function that depends on the connection 
at the finite number of edges of a graph) on A/Q, constructed from a spin network 

for all A & A, which includes distributional besides smooth connections. he-{A) is an element 
of SU (2) associated with the edge and "•" stands for contracting, at each vertex Vk of a, the 
upper indices of the matrices corresponding to all the incoming edges and the lower indices of 
the matrices assigned to the outgoing edges, with all the indices of pk- 

Given a pair a,j the vector space generated by T^j^jj, for all possible contractors associated 
with a, J in the way stated previously, is denoted by T^aAx- Then 

where a, J* run over all the pairs consisting of minimal graphs and irreducible non trivial rep- 
resentation labelings. The sum is orthogonal and the spaces ?iau4 finite dimensional. It 
suffices to define an orthonormal basis within each of them. 

Note that the aforementioned construction of Tiaux holds actually for any diffeomorphism 
covariant theory of connections with compact gauge group. The choice of SU (2) corresponds to 
the case of gravity described in terms of real connection variables. So the generalization we are 
interested in to include both gravity and the electromagnetic field is = Tif^*^^-* ^'^^x^- The 
spin network states for the compound system are denoted by Tajj;p] [A, A], with c, g labeling 
the U{1) coloring and contractors, respectively. 

To extract physical information we will further need a state describing a fiat continuous 
space S at scales much larger than the Planck length, but not necessarily so at distances 
comparable to Planck length itself. States of this kind were introduced under the name of 
weave |^ for pure gravity. Flat weave states \W), having a characteristic length C, were first 
constructed by considering collections of Planck scale circles randomly oriented. If one probes 
distances d » C the continuous fiat classical geometry is regained, while for distances d « C 
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the quantum loop structure of space is manifest. In other words, one expects a behavior of 
the type (iy|gaf)|W^) = Sab + O (^) • It was soon realized that such states could not yield a 
non trivial volume due to the lack of self intersections Couples of circles, intersecting at a 
point, were also considered as specific models of weaves to overcome this defect With the 
recent advances on the kinematical Hilbert space Tiaux it became clear that all proposed weaves 
were afflicted by two undesirable features. First, they are defined to be peaked at a specific 
(flat or curved) metric, but not at a connection. This is in contrast with standard semiclassical 
states in terms of coherent states, for example. Second, the known weave states do not belong 
either to or to a dense subspace of it . It may be possible to come to terms with such 
difflculties by deflning coherent states for diffeomorphism covariant gauge theories or by 
implementing a genuine statistical geometry for instance. Both approaches have recently 
achieved substantial progress. 

Nonetheless, in order to extract physics, there is the alternative possibility of using just the 
main features that semiclassical states should have . Namely, peakedness on both geometry 
and connection together with the property that they yield well deflned expectation values of 
physical operators. An advantage of this alternative is that one may elucidate some physical 
consequences before the fuH fledged semiclassical analysis is settled down. Indeed, such al- 
ternative may be considered as complementary, in the sense of hinting at possible features of 
semiclassical states which could be further elaborated. After its completion, a rigorous semi- 
classical treatment should tell us whether the results arising from this alternative turn out to 
hold or not. The weakness of the treatment resides on its generality, since no detailed features 
of the would be semiclassical states are used -as opposed, say, to the original weave states- and 
hence a set of numerical coefficients cannot be calculated. Evaluating them will be the task of 
the rigorous semiclassical treatment. 

On top of the purely gravitational semiclassical states, a generalization is required to include 
matter fields. For our analysis it will just suffice to exploit the same aspects of peakedness and 
well defined expectation values, extended to include the case of the electromagnetic field. The 
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semiclassical states here considered will describe flat space and a smooth electromagnetic field 
living in it. Such a state is denoted by \W, E_,B_> and has a characteristic length C Since 
no detailed information is used on how the semiclassical state is constructed in terms of, say, 
a graph, as opposed to weave states, the present approach yields results relying only on the 
following assumptions: (i) peakedness of the states, (ii) well defined expectation values and iii) 
existence of a coarse-grained expansion involving ratios of the relevant scales of the problem: 
the Planck length £p, the characteristic length C and the electromagnetic wavelength A. States 
fulfilling such requirements are referred to as would he semiclassical states in the sequel. 



3 The Regularization 

Thiemann has put forward a consistent regularization procedure to define the quantum Hamil- 
tonian constraint of general relativity on Tiaux, both for pure gravity and matter couplings 



31| . The basis of his proposal is the incorporation of the volume operator as a convenient 
regulator, since its action upon spin network states is finite. We use his regularization for the 
Einstein-Maxwell theory, which naturally allows the semiclassical treatment here pursued. 
Consider the electric part of (^. The identity ^ V} = 2sgn{det el)e^^ allows to rewrite 



it as 



^''J 2k^Q^ e^oe^ is ^'^'^^rdetolir.xl- ^ 



Q2e-.oe3 7s 2(detg)3(x) 
^ ^2(detg)i(y)-^^ 



Ihn d'xN{x) ^V{x,e))^E\x) 
X / d'yx.{x,y) {A',{y),^V{y,e)]E\y), (9) 



with Xt{Xi y) = n^=i^(e/2— \x'^ — y"'\) being the characteristic function of a cube with volume 
centered at x and V{x,e) := / d^yXe{x, y)^/detq{y) being the volume of the box as determined 
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by Qab- Remarkably all e dependence resides here. This is possible due to i^Maxweii having density 
weight one and it is achieved at the price of explicitly breaking diffeomorphism covariance. This 
is harmless as far as diffeomorphism covariance is regained once the regulator is removed. This 
is the case indeed |^T|. Next, let S be triangulated into tetrahedra A. Hence, the integral over 
S in (P) is just a sum over the contributions of each tetrahedron A. 

The form of suggests to focus upon the term inside each integral. As we will see below, 
this indeed simplifies the analysis. Let 

©i/] := Jd'xf{x)mx)[Al{x),^V{x,e)^ 

©i/] =: E©k[/]- (10) 

A 

Also let us use the dual ofE_: '■= *E_bc = jj^bcciK^- For a two-surface S, J^e = \ ^bcd dx^A 
dx^ = i JsndK^ ibc dx^ A dx"^ = fsK^nd e, e being the volume two-form. Hence $^(S') := f^e 
is the flux of E_'^ through 5*. Recalling that 

tr (r,K,{^K;,^V{x,e)}^ = tr (r,r„ £ dt r^^'^it) {^^"(.^^(t)), ^V{x, e)}) + . . . 

= £ dt sl'^t) [A-{sl\t)), ^V{x,e)] + ... 

^ -lsUl){A:{sl\0)),^V{x,e)] (11) 

and that, for small tetrahedra, ^^{Fjk) ~ |eafecSj(A)s^(A)£;'', it follows 

f{v) e^^^^^iFjK) tr {r^ {Z^;, V), v'^(.;(A), e)}) 

^ -^/(t;)e^^"6,,,4(A)s^(A)E^l(A) [A^AsT'm. ^V{x,e)] 

= ~mvoliA)E^ [A:,isr\0)),^V{x,e)^ 
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-I / feA{A\x),JV{x,e)\ . (12) 



2 7a 
We have then 

0k[/] = -| f{v) e'''V{FjK)tr (r^ K,^^^ |/.,-\^), ^V{v{A),e)Y) , (13) 

where we are denoting by sj(A), sx(A), sl(A) the edges of the tetrahedra A having v as 
common vertex. As stated, Fjk is a surface parallel to the face determined by sj(A),Si^(A) 
which is transverse to sl(A). 
Hence 

Next one replaces and V{x, e) by its quantum counterparts and adapts the triangulation 
to the graph 7 corresponding to the state acted upon, in such a way that at each vertex v of 
7 and triplet of edges e, e', e" a tetrahedron is defined with basepoint at the vertex f(A) = v 
and segments s/(A), / = 1,2,3, corresponding to s(e), s(e'), s(e") , respectively [^. Here it is 
assumed that e^'^^ eabcSi"'Sj^SK'^ > 0. The arcs connecting the end points of S/(A) and sj(A) 
are denoted a/j(A), so that a loop a/j := Sjoajjosj^ can be formed. Besides, the face spanned 
by the segments S/(A) and sj(A) is called Fjj. 

The action of the regulated operator hereby obtained gets concentrated in the vertices of 
the graph, as expected from the explicit appearance of the volume operator. In successive steps 
we replace 



Qa[N] = Niv{A)) e'^^^^{Fjj,)tr (r^ h,,^^^ [z.-^^^, ^^^(.^(A), e) 



(14) 
(15) 



to obtain 

H'^iN] = — y N(v)(^^] y X 
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X tr I^t' hs^{A) 
xtr It' Kpi^A') 



MNP^E 



MN) 



(16) 



The valence n{v) of the vertex v yields the contribution E{v) = n{v){n{v) — l){n{v) — 2)/3! 
of the adapted triangulation at each vertex of 7. Also, as e — 0, f (A) = t>(A') are the only- 
contributions left over. The final expression for the electric piece of the Hamiltonian constraint 
given in [31| is obtained by the explicit action of this operator on cylindrical functions. We 
refrain from doing that here because the form of the operator (|16D is better suited for our 
approximation given below. 

As for the magnetic part of -f^Maxweii we proceed similarly. Since 



ks 
I) 



1-z / dts''{t)A^{s{t)) + . . . 



I / dta^jK{t)AMt)) + --- 





-i f B^dSa + ... 

-1^''{Fjk) + ... 



and 



f{v)e'''\K^^ - l)tr (t,K,^a) \A^} 



t-volisj, SK, sl) 5: fiv) B\v) { A\{v), V^(x, e] 
.3! 



i-vomf{v)B\v){A\{v),^V{x,e] 
.3! 



I- fix) B\x) A \ A\x), ^V{x,e) 



(17) 



;i8) 
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we can write 



where 



(19) 



(20) 



Contrasting S}^ with 0}^ we notice (i) different numerical factors: for the former 
while — 1| for the latter and (ii) to leading order, the magnetic flux €"^^^{11^ — 1) has as its 



counterpart the electric flux e'^^^E, 



OiJK ■ 



The quantum version of the above operators is obtained using 
Finally, one gets the regularized magnetic piece of the Hamiltonian constraint as |^| 

^7p(A')'V^) (^A/iv(A') 



X 

u(A)=t.(A')=?^ 

lU 



X e 



4 il'Sp(A') 



1 . 



(21) 



The electric and magnetic pieces of -f^Maxweii can be treated in a unified manner in terms of 
fluxes. To see this recall that for abelian gauge fields 



where ^^{Fjk) is the flux of the magnetic field through the surface Fjk- 
Then the full electromagnetic Hamiltonian is 

1 / V N \ 

-f^Maxwcll[^] = + 



X 
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^ JKLMNP 
Ac c 



tr { Tjh 



i'l'Sp(A') 



1 - <|.^(Fj^)|.^(Fl 



MN) 



m) 



Let us emphasize the structure of the above regularized Hamiltonian. There is a common 
gravitational factor included in the SU (2) trace. The basic entities that regularize the electro- 
magnetic part are the corresponding fluxes: one is associated with the magnetic field, which 
enters through a product of exponential flux factors, while the other is related to the electric 
field, entering in a bilinear product of fluxes. Thus, in the quadratic field approximation the 
effective Hamiltonian preserves duality invariance. Nevertheless, the magnetic sector includes 
higher powers in the field expansion. Hence, nonlinearities in the field equations, inducing 
duality violations, arise only via the magnetic field. 

Before proceeding let us recall that, according to Thiemann's conventions, in flat space we 
must have 



Maxwell 



1 



2Q^ 



(23) 



where Q is the electromagnetic coupling constant. The electromagnetic potential is denoted 
by and the corresponding electromagnetic tensor by F^^. The units are such that the 
gravitational connection has dimensions of 1/L (inverse Length) and the Newton's constant 
hi has dimensions of L/M (Length over Mass). Also we have that [\E\/Q'^] = M/L^. Taking the 
dimensions of to be 1/L, according to the corresponding normalization of the Wilson loop, 
we conclude that [E] = [B] = and [Q^] = 1/{M L). In our case we also have [h] = M L, 
which in fact leads to aEM = Q^Ti to be the dimensionless fine-structure constant, as defined 



by Thiemann ||3TI] 
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4 General structure of the calculation 



The effective Maxwell Hamiltonian is defined by considering the expectation value of the U{1) 
gauge sector of the quantum Hamiltonian constraint with respect to \W,E_,B_). Inside this 
expectation value operators are expanded around all relevant vertices of the triangulation in 
powers of the segments s1{A), which have lengths of order ip. In this way, a systematic 
approximation is given involving the scales ip <K £ < A, where A is the De Broglie wavelength 
of the photon. Our corrections to the Maxwell Hamiltonian arise from such an approximation. 

We do the full calculation of the magnetic sector, including the non-linear contributions to 
order £p. Next, to obtain the electric sector, it is enough to consider only the quadratic terms 
in the magnetic Hamiltonian and make the replacement B ^ E. 

In the case of the magnetic sector, the general form of the expectation value is (recaUing 
that = n^/^p ) 

(^ai . . . F^^{v)) t,...a^P'^P''''-''"'-{v, s{A), s{A'))\W, E,B>. (24) 

To proceed with the approximation we think of space as made up of boxes, each centered 
at a given point x and with volume ~ x. Each box contains a large number of vertices 
of the semiclassical state {£■ » ip), but is considered as infinitesimal in the scale where the 
space can be regarded as continuous. Also, we assume that the magnetic operators are slowly 
varying inside the box {ip « A^^), in such a way that for all the vertices inside the box one 
can write 

{W, E, R\ . ..F^,{v) ...\W, E, B) = iiea,cB%x). (25) 

Here — daA^ — dbA^, R"'{x) is the classical magnetic field at the center of the box and is 
a dimensionless constant to be determined in such a way that we recover the standard classical 
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result ( p3D in the zeroth order approximation. In the next section we show that 

/^=hr ' (26) 



'-p. 

with T being a parameter defining the leading order contribution of the gravitational connection 
to the expectation value. Applying to (|^) the procedure just described leads to 

2 



Box jjgBox 



3! EM 



1 - 1 



E V < ^' ^' 5|-^r„,....„^''^^'-P" (A), s{^'))\W, E,B> 



pqpi q\...pn Qn 
•dm 



Box ^ ^ 

= / d'x iV(x)^,,,,(f) . . . (^"^ . . T,,...,„^'^^'^''-^'"'?"(x). (27) 

The box-averaged tensor T(ij...a^'"'^i''i-"P"(x), defined by 

rp p5pigi...p„g„/-jN ^ V 2 8 \ ^ / 1 

^ ^ ..Voxl3!i5(.); 1 

X < W^, I, S|^t,...,„^^^^''-P"''"(t;,s(A), s(A'))|l^, E, B >, (28) 

is constructed from fiat space tensors like 6ab, ^abc- In this way we are demanding covariance 
under rotations at the scale C. 

When averaging inside each box, the scaling of the expectation values of the gravitational 
operators is estimated according to 

{W,E,R\...Aa...\W,gB)^...U^-^] (29) 



{W, E, B\ . . . ^V^ . . .\W, E, B) ^ . . . . . , (30) 

respectively. In our previous work |^ we have set T = on the basis that the coarse graining 
approximation, defined by the scale C, does not allow for the connection to be probed below ^. 



On the other hand, by adopting naive kinematical coherent states for representing semiclassical 
states, one would set T = 1 for two reasons: (i) to guarantee that ( PUD yields just zero in 
the limit ^ — , in agreement with a flat connection and (ii) because such an scaling would 
saturate the Heisenberg uncertainty relation Nonetheless, physical semiclassical states 

may imply a leading order contribution with T 7^ 0, 1, thus we choose to consider T as a free 
parameter here. 

In order to make the transition to the electric sector it is convenient to express the effective 
Hamiltonian in terms of the magnetic field, which amounts to a redefinition of the expression 
p7| ) in the form 

H^ = y' d'xN{x)B,^{x)...B,Sx) i?a,...a„''^-""(x). (31) 

The relation between the box-averaged tensors R and T is 

rri pqp\qi...pnqn(z;\ _ ^pq ,pngn p rri...rn(z;\ /ooX 

■^ai...am K-^J ~t r-j . . . t r„ -"'ai...am \-^J- 

By expanding ( ^21) to different powers in Sj(A) one can systematically determine all possible 
contributions to the effective electromagnetic Hamiltonian at a given order in ip. 



5 The Calculation 

In this section we provide the details of the calculation of the Maxwell effective Hamiltonian 
up to order £p. Let us start with the magnetic sector. 

The two main ingredients in (^2]) which contribute to the expansion in powers of the seg- 
ments s/(A) are: (i) the trace factors involving the gravitational operators and (ii) the magnetic 
flux through each surface F/j(A). 

First we calculate the flux of the magnetic field through the surface Fjj. A convenient way 
to do this is via the Stokes theorem 

$^(F7j) = / Ban'' (fx = f fv X 1) d^x 
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dt s''{t)Aait) 

aij 

v+Sj rv+Sj pv 



nv+Sj rv+s J rv 

= / Aadx''+ ' Aadx" + I ^ A^dx"". (33) 

Jv J v+sj J V+Sj 

Here the notation is s} = (s"} and analogously for v. We are using straight line trajectories 
joining the vertices of the corresponding triangle. 
The basic building block in (|33|) is 

Mx)dx'' = / M^l + t {V2 - Vi)) {V2 - ViY dt 



vi JQ 



Mvi + t A) A" dt 

Jo 



1 + + l^iA%f + . . .) A^A^iv), (34) 

with A*^ = {v2 — Vi)"'. The infinite series in parenthesis is 

F{x) = 1 + + + + • • • = (35) 

yielding 

MS) dx^ = F{A^da) (A^Mvi)) . (36) 

In the following we employ the notation A" Va = A ■ V. Using the above result in the three 
integrals appearing in (^) and after some algebra, we obtain 

^""{Fjj) = F,{si-V,sj-V)s'}sUdaA,iv)-dtM^)) 

= F,{sj-V,sj-V) s"js'jeabcB'{v), (37) 

where the gradient acts upon the coordinates of v. The function Fi is 

xy[y-x) ^^[n + l)\ x-y 

Let us emphasize that Fi{x^y) is just a power series in the variables x and y. Expanding to 
fourth order in the segments we obtain 

^^{Fu) = (l +l{s'i + s'j)d,+^{s'jsj + sjsj + s'ysj)d,d, + ..}jx 
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x]^s'\s%,,B\v). (39) 



Notice that the combination 



\s'\s)eabc = An^ (40) 



is just the oriented area of the triangle with vertex v and sides Sj, s}, joining at this vertex, 
having value A and unit normal vector ric- 

In order to make the bookkeeping clear, let us denote by T the combination arising from 



([2^ ) whose expectation value we are calculating 
1 

2Q 

Here 



T = --^^mL^mPA' e^^^ 6^^^ (e-4^(^-) - l) (e-*^(n..) _ i) . (41) 



(42) 



WiLA = tr[Tih,^^A) 



\ V 



Some remarks are in order before we proceed further. Our final goal is to obtain a power 
expansion of ( |H| ) up to order £p. Since, as we will show in the sequel, the normalization factor 
converting magnetic operators inside the semiclassical expectation value into classical fields 
outside the expectation value is proportional to (£p)~^, we have to take some care regarding 
the expansion of the given quantities in powers of s". A detailed power counting analysis in 
the expression ( p]) shows that the term Wi^A) in (0) is to be expanded up to order s^, while 



each magnetic factor (^e~^^^ ^^"^^ — 1^ is required to have the following properties: the terms 
proportional to F_ are to be expanded up to order s^, those proportional to up to order and 
finally those proportional to up to order s^. This will lead to the following contributions in 
T: the terms proportional to include the expansion up to order s^,the terms proportional to 
include the expansion up to order and the terms proportional to include the expansion 
up to order s^^. The final result is that the semiclassical expectation value of the magnetic 
contribution T will be proportional to £p, which is incorporated in the volume element, times 
corrections up to order £p. 
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Let us now continue with the calculation of the contribution to (^) due to the magnetic 
flux by writing the expansion 



n=l ^• 

+0 (43) 



where 



MijK{A) 








M2JK(A) 






MsjK(A) 


■= S'ks^j^ix] + XjXK + 

1 




Fat - S'^K 




+ ^^^^EabiXK + Xj) 


Fed 


— a"- Q 



-cdy 



{Xj + XK)FabFcd 



4-3! 

i 



^ _ gj^afeZcdZe/; (44) 

according to the previous analysis. We are using the notation x/ = s/- V = da- Let us remark 
that, contrary to the electric case, the magnetic contribution will incorporate non-linear terms 
due to the expansion of the exponential in powers of B. This implies that the exact duality 
symmetry of Maxwell equations in vacuum will be lost due to quantum gravity corrections. 

Next let us consider the gravitational contributions to (p2D, arising from (^2]), which we 
expand as 

WiLA = slWia + Sls'lw^ab + sls'lslWiabc + 0{s^w), (45) 

with 



w. 



-- ^[Aia, W], Wiab = ^eijk[Aja, [Akb, W]], Wiabc = -^^i^ja, [Ajb, [Ac, V^]]]- (46) 

The scaling properties of the above gravitational operators under the semiclassical expectation 
value are 

f 3/2 /f X""^ 

{WEB\ . ..w,a,...a^ ■ . ■ \WEB) ^ ^ ( f 1 • (47) 
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Let us emphasize that the result p7| ) is a consequence of the scahng of the expectation value 
of the connection given by 

For the product WnAWiPA' we need only 



with 



U2LP 

UsLP 
UiLP 



WiLAWiPA' = U2LP + U^LP + UiLP + O [s^tV^^ 



sls'pWiaWid, 

S i^S pS pWiaWide + S j^S j^S pWiabUlidi 

"" ~\~ Si^Si^SpSpWiab'UJiiie~\~ Sj^Sj^Sj^SpWiabc^id- 



(4J 



S j^S pS pS pWiaWidef 

Here all the w's are evaluated at a common vertex v. 

At this level it is convenient to state the result (no sum over L) 

si si Wiab = I si si eijk [Aja, [A^b, W]] 
o 

= g ^2 4 ^ijk {AjaAkbW - AjaVVAkb - AkbVVAj 
= 0, 

which holds due to symmetry properties. This leads to 

UsLP = 0. 

After taking the expectation value the terms contributing to order £p in (p]) read 



(49) 



VAkbA.a 



^ja 1 



T = To + Ti+T2 + 0(^4) 



Tn 



1 ^ ^JKL ^MNP 



2Q2 £4 



-'- -'- ^JKL MNPtt 

'W^P 

1 1 JKL MNP 



2LP 



U2LP MijK M[^^ 

MijK M^MN + M2JK M[m^ 

U2LP {MuK M^A^JV + MsjK M[mn + M2JK 
+U4LP MijK M[j^f^ 



2MNj 



(50) 

(51) 

(52) 
(53) 

(54) 
(55) 
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Now we are ready to calculate the different contributions to the magnetic sector of the 
Hamiltonian (E^), which we parameterize in terms of the tensor Wa^,,,aJ^'^^"'^" introduced in 



Recalling that we are only interested in the pieces which are symmetric in the indices 
ri, r2, . . . , r„, the contribution Tq produces 

2 ,2 



p nra _ ^ [ ^ 8 ] JKL MNP 



X 



X s'i.s'jsl s'l,s'ls'p {WEB\w^cWid\WEB). (56) 

In order to simplify the product of vectors (s'm^) appearing in the sequel and also to exhibit 
the internal symmetry properties of the quantities involved, it is convenient to keep in mind 
the relations 

e^^^ s^j si = det{s) e"^^ det{s) = det{s\), e'^^^e^fe, = 25^ (57) 
In this way, ( ^6|) can be rewriten in the simpler form 

D rira ^ V- _J_ ( ^ ^ ^ //^ 

° 2g2 \^\E{v) fp 

rig Box ^ \ V // v(A)=v(A')=v ^ 

X det{s) det{s') {W, E, B\-{wi^\ wp} \W, E, B). (58) 



The above equation implies 



1,2 P / P \'^^ 1 

- 2Q^£^^^£^l,Zj -2g^^ ' ^^^^ 



which reproduces the zeroth-order magnetic contribution (p3D with the choice 
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Now let us consider the correction arising from Ti, which leads to the following contribution 
in the effective Hamiltonian 



Box(a;) v€ Box \ V 



v{A)=v{A')=v 



In the above expression we have interchanged the summations over A and A' in order to 
rewrite WicWid as |{wic , Wid}- We further separate the above contribution to Hf^^ in two parts: 
(i) the first contains two powers in the magnetic field and leads to Rnai^^^ and (ii) the second 
one contains three powers in the magnetic field and leads to the completely symmetric tensor 
j^^_^rir2r:i _ gjnce we have no symmetric tensor with three indices at our disposal, the latter 
contribution is zero. Thus we concentrate in the first one 

Srs'j s'ls'l, s'Z e'\b e\g{W, E, B\{wic , w,d} \W, E, B). (62) 

In this case the internal symmetry properties are hard to make explicit. In order to determine 
whether or not the above contribution is zero we contract (p2| ) with the only three index tensor 
at our disposal: The result is 



x{W,E,B\{w,,,w,d}\W,E,B) [4s'Ijs'Is\m + 4s'%s'l,s' 



qM 



(63) 



Inside the expectation value, upper spatial indices have been lowered by the flat metric. By 
symmetry requirements, the second term in parenthesis in (K3) yields zero. Nevertheless, the 
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first one gives the result 

p aim _ ^ ^ p7 ( ^airri ff-A^^ 

- ''^Q^^^Z^ ' ' ^^^^ 

which produces a parity-violating term in the magnetic sector of the effective Hamiltonian. 

The next contribution arises from T2 and can be separated into three pieces 

/ \ 2 

^ I ^ \ 0^ ST^ ^ ^JKL MNP 



- -;r e e X 



X {W, g B\ U2LP {MuK M'^j^^ + M^jK M[j^^) \W,E,B), (65) 

x{W, EB\ U2lpM2jkM'^^^ \W, E, B), (66) 

X {W, E, B\ U^LP MuK M[j^^ I W, E, B) . (67) 
Let us start discussing Hf^^. After some algebra we obtain 

ttB _ V — V _J^^JKLMNP 



1 1 



-55,45^5>^5}^F,,Z„,F,,) F,, \W, E, B), (68) 
which naturally splits into the following pieces 

ttB _ \^ 1 f ^ \ p3 '^_JKL^MNP^ 
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X sl4 s'nS% {W, E, B\ WiaWid (s'ks'j^{x^j + xjXk + x\)F,}j F^g \W, E, B), 

(69) 

ttB _ 1 f ^ \ p3 \^ j_ JKL MNP ^ 



X sl4s%sl {W, E, B\ {s'ks'jsIs^j 



-{xj + XK)F^tF^ 



+j^F^tixK + xj)F^ 



4-3! 

F,^\W,E,B), (70) 



1 /2 8 \\, ^ i 



X sl4s%sl {W, EB\ WM [s^js'ks^sls^j^F^.F^X^,^ F^^ \W, EB). 

(71) 



Next we calculate the corresponding R tensors. The first one is 



2 - ,2 



eNKj (s'^) e-'^^ si {2sJ sj + sj s^^) {W, E, B\ w,aw['\W, E, B), 

that we have constructed using the most general tensor with four indices, which is symmetrical 
in oi and 02- The terms proportional to contribute with B^d^-d^Bi, to the effective Hamil- 
tonian. Keeping this term will be useful to construct the electric sector later. Notice that we 
have also introduced the inverse matrix {s~^)^ such that 

{'-fm'Q-^Q^ 'm{s-')^-SI (73) 

which elements scale as ip^. 
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Our next contribution, arising from (|70D, is zero because the object (no sum over J and K) 

SkSjS^s} Z„„ {xj + XK)F,t, (74) 

which is symmetrical in J K, appears contracted with e'^^^. 
Now we consider 

2 4 



Box(x) t>eBox v{A)=v{A')=v 



X 

y<{W,E,B\w^aW,^\W,E,B), (75) 
where we have introduced the relation 

emfec s J = rfet(s) exjL (s^^) ^ , (76) 
obtained from (|57D. The further properties 

e^^^ euKj^VKj^xKj =^ L = U = V = X, e^^'^ ^ikj^ikj^ikj = 2e^^^ ei23ei23ei23 = ~2, 

(77) 

allow us to rewrite the corresponding tensor in the simpler form 

= 5 W {1^.^) ^''^''^ ^^''^'^ 

tJSBox ^ \ V / / v(A)=v(A')=v -r 

^ ('"')ai ('"')a, ('"Oaa ^' ^1 1^' ^' ^'^^^ 

This shows explicitly the symmetry in the indices ai, 02, as and leads to 

-^21801020304 Q2 £7 ~£2 ^"^^ '^^ '^0102 '^0304 ~l" '^0103 '^0204 ~l" '^0104 '^0203^ ; ('''9) 

where we have written the most general four index tensor at our disposal, which is completely 
symmetric in three indices. This implies the complete symmetry in the four indices. The 
correction (|79D leads to cubic non-linear terms in the equation of motion. 
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Now we continue with the correction arising from (^), which reduces to 

TB_S-'^f'^^\/]3 4 ^JKL MNP 



\2 



?;GBox ^ \ V / / v(A)=v(A')=v \ / P 



e e X 



2( 

xsjs^ (ly, £ 5| w^aWid (sls^jXjF^y) {s%s%x'j,jF,,) \W, E, B), (80) 



because the terms quadratic in F_ are symmetric in J, K and M, N respectively. 
/^From the above we read 



2 ^ ,2 



V3!i^(^); .(A)irA')=^ (3!)2£^ 
y<{w,E, B\ wm\W,E, B). 



Here the internal symmetry properties are rather obscure. Nevertheless, the symmetry 
properties induced by the classical magnetic field factors imply that the above tensor must be 
proportional to the most general tensor with four indices which is symmetric in ai 02- Then we 
have 

iP' P ( ( \'^'^ f \ 

This contribution is of the same kind as the one given by i?2ii"^"^'^^^- 
Finally we are left with 



Box(:f) DeBox(x) ^ V-'--^V'^;/ v(A)=v{A')=v 

V ^JKL MNP u V Jx ly f Jajbjc'^ 



2Q2 \3\E{v) ^ 4. 

^ \ y ' / v(A)=v(A')=v 

[si44s'^) e,/^ e,/(H^, E, B\ {w,d, w,abc} \W, E, B), 

(83) 

which leads to 



2 ,2 



i)GBox(^) ^ \ V^;/ v{A)=v{A')=v 



X (s'-i)' (s'^ 4 s'^ ) (l^, I, B\ {wr, Wiabc} \W, E, B). (84) 
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Taking the symmetric part, we have 



-K23 



^5) 



Adding all previous contributions, we obtain the magnetic sector of the effective Hamilto- 
nian, up to order £p. 



H 



B 



2T 



B + 



^6) 



The numbers 9i are linear combinations of the corresponding /t^ appearing in the tensors R. 
The correspondences are 



l^lt I^U) ^2, ^3; l^i Kii 6j, , Ks —>■ 9s- 



As it is pointed out after (P^, the electric sector of the effective Hamiltonian can be obtained 
by changing B into E in the quadratic contribution of the corresponding magnetic sector. In 
this way the complete electromagnetic effective Hamiltonian becomes 

/ n \ 2+2TN 



H 



EM 



+02 tpE^dadbW + Os^P [ B-{V X B) + E-{V X E) 



2 n2 



2T 



B^ + 



up to order ip. 
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6 Modified Maxwell equations and dispersion relations 



in vacuum 



Since no confusion arises in the sequel, we eliminate the underline in all electromagnetic quan- 
tities. From the effective Hamiltonian (^) we obtain the equations of motion 

A(V xB) — + 24^3 V2(V xB)- 2esipV^B + AO^C^ — ilV x (B^B) = 0, (89) 

c ot \lpj 

1 r)R 

A{V xE) + -— + 24^3 V2(V xE)- 2esipV^E = 0, (90) 
c ot 

where 



2+2T 



(91) 



The above equations are supplemented by the condition V ■ -B = 0, together with the constraint 
V ■ -E = 0, appropriate for vacuum. 



Next we calculate the dispersion relations arising from the modified Maxwell equations ( [89| ) 
and (pOl). Neglecting the non-linear part and introducing the plane wave ansatz 



(92) 



we get 



Eo-k = 0, Bo-k = 0, 



UJ 



(kxEo) A-2e:i{£pky -2iesip k'Eo--Bo = 



(kxBo) A-2e^Up kV -2i^8£p r^o + -Eq = 0, 

L -I C 



(93) 
(94) 

(95) 



which imply the following dispersion relations 

2+2T 



u = ck 



2 03 {kipf ± 208 (kip) 



(96) 
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The ± signs correspond to the different polarizations of the photon. From the above we obtain 
the speed of the photon 



V 1 duj 
c c dk 



= 1 ± 4^8 (Hp) - Qe^{Mpf + e^ {Mp)™ + ... (97) 

C=l/k 



The last expression gives v expanded to leading order in £p where we have estimated C as 
1/k, which is its maximal value. Clearly (|97D is valid only for momenta satisfying (ipk) « 1. 

To first order in (kip) we have only the helicity dependent correction found already by 
Gambini and Pullin 0. As far as the T dependent terms we have either a quadratic (T = 0) 
or a quartic (T = 1) correction. The only possibility to have a first order helicity independent 
correction amounts to set T = —1/2 which corresponds to that of Ellis et. al. ||10[- However, 
we do not have an interpretation for such a value of T. 

7 Effect of non-linear terms 

In this section we explore some implications of the non-linear term in the Maxwell equations 
induced by the quantum gravity corrections to order £p. 



Following reference |^ we study the propagation of waves in the presence of a constant 



magnetic field Bq. To do this, let us write 

B = Bo + b (98) 

and consider only the contribution of the non-linear term in the effective Hamiltonian. Namely, 
let us set 03 = 9t = 9s = 0. After hnearizing in b, the field equations reduce to 

1 dE 

V xb--— + 46^ C^fpV x{B^b + 2{Bo -1^)130) = 0, 

V-E = 0, V-6 = 0, Vx^ + -— = 0, (99) 

c ot 
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with 9a 



2T 



6*4. Now we look for plane wave solutions of ( pQ]) with 



(100) 



obtaining the conditions 



UJ- 



Eo-k = 0, bo-k = 0, {kx Eo) 60 = 



UJ 



1 + 404 C'tp B^) {k X 60) + -^0 + 8 04 C'tp {bo ■ Bo) {kxBo)=0. 



(101) 
(102) 



Substituting in ( |102| ) the expression for 60 obtained from the third equation in ( |101|) we are left 
with 



- A;2 (1 + 4 ^4 CHl Bl) \eo = -S 0^ CHl I {k x Eq) ■ Bo)) (k x Bq). 



(103) 



Since Bq and k determine a plane, it is natural to study separately the propagation of waves 
with polarization parallel and perpendicular to this plane. We will express the answer in terms 
of the refraction index 



n 



k c 



UJ 



(104) 



/^From ( |103| ) we obtain the following refraction indices {h = 1 = c) 

^ 1 



77.11 



1 - 204 



2 n2 
2 -"O' 



UJ 

1-204^ fpB^ (1 + 2 sm\ 
ui'^ 



(105) 



for parallel and perpendicular photon polarization, respectively. Here is the angle between k 
and Bo- 

These results can be compared with similar effects in quantum electrodynamics p9 



= 1 + 0\\Blmi^(l), 
n_L = 1 + 0^Blsm^(j). 



(106) 
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where 

2e^ 7e^ 

As expected, quantum gravity induced effects are much smaller than purely quantum elec- 
trodynamics effects, but the former present different signatures. In particular, the indices 
arising from quantum gravity are frequency dependent and also nii is independent of (f). 



8 Photon time delay 



Notice that our considerations assumed a coarse-grained fiat spacetime rather than a FLRW 
model. The latter would seem more appropriate for GRB traveling cosmological distances. 
In particular, it would be interesting to study the redshift effects in the photon time delays 
induced by the energy dependent corrections to the velocity. We are able to estimate these 
effects as follows. 

Using the flat FLRW metric let us calculate the present time delay of two photons emitted 
simultaneously with different momenta and hence different velocities. We set r — — t — ti 
as the emission coordinates and r — 0,t — to as the detection coordinates in the comoving 
cosmological system t,r,9{t) — 0. The definition of the velocity V{t) — R{t) with ^ < 0, 
leads to 

rr-i fto V(t) 

'■' = /. ^"-"i. sliy*- 

Here R{t) is the scale factor in the FLRW metric. Notice that the above equation leads to the 
standard red shift result for photons moving with V — c. On the other hand, quantum gravity 
corrections predict 

-^1 + 9 (ipk) (109) 
c 

to leading order. We are interested in discussing the situation where two photons are emitted 
at ri with different velocities Vi and Vi + 5Vi, and arrive at r = at times to and to + 5to 
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respectively. Then we have 

Subtracting the two expressions for the fixed coordinate ri, we obtain 

In the fiat FLRW universe we have 

_ to _21,_1 1, R{t) _ 1 

'-(TTi)^' '°-3^' ^'--^(TT^^"' ^-IT^' ^^^^^ 

where 2; is the red shift. This leads to 

^to = ^ r ..f\,,, SV{z), (113) 
cHq Jo [1 + zy/"^ 

where zi is the red shift of the source and 6 V{z) is the difference between the velocities of 
the two photons, arising from the quantum gravity corrections ( |109|) . Using the zeroth-order 
relation for the red shift of the photon momentum, 

ki.) = ^*„. (114) 

we obtain 

= mo ip = 9il + z) i5k)o ip, (115) 

where {6k)Q is the present day observed momentum difference between the two photons. 
Substitution of the above equation in (|113|) yields 

6to=(^] moip Ui+z^y/'-i). (116) 







The above result differs from the corresponding one obtained in the second and third papers 
ofRef. 0. 
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9 Summary and Discussion 



In this paper we have considered the propagation of photons in a semiclassical background 
provided by loop quantum gravity. An effective electromagnetic Hamiltonian, given by (^), was 
identified with the expectation value of the electromagnetic piece of the Hamiltonian constraint 
for the Einstein-Maxwell theory with respect to a would be semiclassical state. The state used 
was assumed to approximate a classical flat metric, a classical flat gravitational connection and 
a generic classical electromagnetic field, at scales larger than the coarse-grained characteristic 
length C, where C » £p (the Planck length). To leading non trivial order in ip, photons of 
wavelength A, where £ < A, acquire Planck scale modified dispersion relations as compared 
to those in classical fiat spacetime. This in turn yields the effective speed of light (0) which 
involves two types of corrections. One of them is just that of Gambini and PuUin including 
the helicity of the photon, whereas the other depends on the scale C Moreover, the latter 
type contains a parameter T that encodes the scaling of the gravitational connection under the 
semiclassical expectation value. 

When estimating the coarse-grained characteristic length hj C = 1/k, which is its maximal 
value, the following values of T are prominent: (i) T = can be understood as that the 
connection can not be probed below the coarse graining scale C The corresponding correction 
scales as (kip)'^. (ii) T = 1 may be interpreted as the analog of the coherent states analysis, 
where such states saturate the Heisenberg uncertainty relation inside a box of volume C^: 
Ag ~ Av4 ~ ^ and AqAA ~ # Hf. Then the correction behaves like (kip)^. 



m 



Interestingly, a value T = — | leads to a helicity independent first order correction (i.e. (kip)) 



similar to that of Ellis et al. |]T0|- We do not have an interpretation of this case though. 

A prime candidate for testing the effects which are linear in the energy would be the Gamma 
Ray Bursts that travel cosmological distances and which might be detected with a time reso- 
lution beyond 10~^ seconds. This seems possible in future spatial experiments 

Moreover, new non-linear terms in the Maxwell equations appear. These terms are not 
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present either in |^ or |T0|. We have explored the significance of this contribution to the 
propagation of photons in a constant strong magnetic field. The corrections obtained in the 
corresponding refraction indices are much smaller than similar effects in Quantum Electrody- 
namics. Nevertheless, quantum gravity corrections have distinct signatures: a main difference 
is that the speed of photons with polarization parallel to the plane formed by the background 
magnetic field and the direction of the wave is isotropic. 

Our results should be taken as first steps in the exploration of possible observable conse- 
quences of quantum gravity. We have given evidence that dispersion relations of the form (|l|) 
can have origin in the microstructure of spacetime. It is expected that the recently proposed 
coherent states for quantum gravity and gauge theories in and/or statistical geometry 



37|, will come to terms systematically with the unknown numerical coefficients we have left 



undetermined in our calculation (See also [0].) Interestingly, a quantum field theory in effec- 
tive spacetimes might possibly emerge along these lines. Yet another avenue in the context 
of canonical quantization of gravity and gauge theories has recently emerged to understand 
the semiclassical regime. It is aimed at establishing a relation between Fock space and TYaux 

Further work remains to be done in the framework here developed. For instance, in the 
case of inflationary cosmology as well as in the study of the Hawking effect, use is made of 
scalar fields with non standard near Planckian frequency dispersion relations to model the 
effect of short distance physics on the quantum fields |^. Indeed a systematic study of the 
modifications induced by quantum gravity along the lines we have developed could be performed 
to investigate whether the dispersion relations used in can be accounted for. 

Finally, we stress that the dispersion relations we have found, as well as those in [p, |10|, ^ 
pl| are Lorentz symmetry violating. This is not necessarily an issue as it has been extensively 
discussed previously [||, In fact they may alleviate some long standing astrophysical and 



cosmological problems Remarkably, there has been considerable progress in setting 

bounds to Lorentz invariance violation [^, ^ and to the values of some coefficients in the 
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effective Maxwell equations 
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